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Abstract 

Relativistic formulas for the electric quadrupole moment of a hydrogenlike atom, induced 
by the hyperfine interaction, are derived for ns and npi/2 states. Both the magnetic dipole 
and electric quadrupole hyperfine interactions are taken into account. The formulas are 
valid for ions with arbitrary nuclear charge and spin. The induced quadrupole moment is 
compared with the nuclear quadrupole moment for a wide range of hydrogenlike ions. 

1 Introduction 

In Ref. [1] it was first noted that a hydrogen atom in the triplet s state has an electric 
quadrupole moment. This is caused by the admixture of the d 3 / 2 states to the s state due 
to the hyperfine interaction between the electron and the magnetic dipole moment of the 
nucleus (HFI M ). In Refs. [2-4], the quadrupole moment (QM^) of a hydrogenlike atom in 
the Is state was calculated in the full relativistic theory: 

Q M (l Sj ^,I, / u) = ^±ll_^Q 1) (!) 

Here 

Qi = \c?^a\ (2) 

Z is the nuclear charge number, / is the nuclear spin, a is the fine structure constant, a is the 
Bohr radius, m e and m p are the electron and proton masses, respectively, \x is the magnetic 
moment of the nucleus, /ijv = \e\H/(2m p c) is the nuclear magneton, 7 = [1 — (oZ) 2 ] 1 / 2 , 

_ = , 1)F -i-m F(2F-1)(2F + 1) 

1 { ' /(/ + F + l/2)(/ + F + 3/2) 1 ' 

and F is the total atomic angular momentum. For the ns state with F = 1, I = 1/2 the 
corresponding calculation to the lowest order in aZ yields [5]: 

Q M (n 8 ,l,l/2,/x) = -^-^Qi. (4) 

A hydrogenlike atom has the QM^ also in p\/2 states, which is caused by the admixture 
of P3/2 states. The relativistic theory of the QM M for a hydrogenlike atom in the 2pi/ 2 state 
with F = 1, I = 1/2 gives [6]: 

Q„(2 Pl/2 , Z, 1/2, M ) = - J_^i(4iV 3 + AT 2 + 9iV - 6)^Qi, (5) 
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where N = (2 + 27) x / 2 . 

The electric quadrupole hyperfine interaction of the electron with the nucleus (HFIq) can 
also induce a nonzero electric quadrupole moment of the atom (QMq), mixing the s states 
with the d states, and the pi/ 2 states with the p 3 / 2 and the f 5 / 2 states. However, since for 
light atoms the HFIq is much weaker than the HFI^ and, in addition, it can contribute only 
for atoms with the nuclear spin / > 1/2, it was generally disregarded in previous calculations. 

In the present paper, we derive relativistic expressions for the electric quadrupole moment 
of a hydrogenlike ion being in ns and np\j 2 states. The derived formulas are valid for ions 
with arbitrary nuclear charge, spin, magnetic dipole and electric quadrupole moments. Both 
magnetic dipole and electric quadrupole hyperfine interactions are taken into account. 

Relativistic and Heaviside charge units (h = c = 1 , a = e 2 /An) are used in the paper, 
the charge of the electron is taken to be e < 0. 



2 Induced quadrupole moment 

The hyperfine interaction operator is given by the sum 

Hum = + Hq, (6) 

where and Hq are the magnetic-dipole and electric-quadrupole hyperfine interaction 
operators, respectively. In the point-dipole approximation, 



|e| (o-Jgx r\) 

H » = ^ 3 > (7) 



and, in the point-quadrupole approximation, 

m=2 



H Q = -a Q2njf 2m {n) . (8) 
m=-2 

Here the vector a incorporates the Dirac a matrices, \x is the nuclear magnetic moment 
operator acting in the space of nuclear variables, Q2m = J2i=i r i^m(fii) is the operator of 
the electric quadrupole moment of the nucleus, i] 2m = C 2m (n)/r 3 is an operator that acts on 
electron variables, ft = r/r, r is the position vector of the electron, fj is the position vector of 
the ith proton in the nucleus, Ci m = \jA.%j (21 + 1) Yi m , and Yj m is a spherical harmonic. It 
must be stressed that the electric quadrupole interaction should be taken into account only 
for ions with I > 1/2. Thus the quadrupole moment of a hydrogenlike atom in the \ A) state, 
induced by the hyperfine interaction, is given by 

E V^ A (A\Q zz \N)(N\H u + Hq\A) , s 
Q,(A) + Q Q (A) = 2 E A -E N ' {Mfa = Fa) ' (9) 

where Q zz = —r 2 (2>n 2 z — 1), \A) and \N) are the state vectors of the total (electron plus 
nucleus) atomic system, Ea and En are the related energies. 

Let us consider first the ns state. Integrating over angular variables in © yields 

Q^ns, Z, I, y) = Y^w(n - l|r 2 |£i(2; n, -1)), (10) 

and 

Q Q (ns, Z, I, Q N ) = aTQ N (2{n - l|r 2 |£ 2 (2; n, -1)} + 3<n - l|r 2 |£ 2 (-3; n, -1)}). (11) 
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Here (i = (II\/j, z \II) is the nuclear magnetic moment, Qn = 2 (II IQ20 1 H) is the electric 
quadrupole moment of the nucleus, 

£ , F = / + l/2,//l/2 

W , F = I -1/2,1*1/2 , (12) 
,7=1/2 



r = < 







|£i(ft';n,/«)} = 



^2(/e';n,K)) = X] 



n' n')(n' n'\a x r 2 \tik) 

En,K E n > ,k' 

n' k'} (n' n'\r~ 3 \nK) 
~E — E 



n,K J-'ri'/K' 



(13) 



(14) 



o x is the Pauli matrix, the vector |ra/e) = ( r ^, raK J ((n«|n/«) = j^(g 2 riK + f 2 K )r 2 dr = 1) 

V fJriK ' 



consists of the upper and lower radial components of the Dirac wave function defined by 



\nnm) = 



gnK,(r)9, Km (n) 

'/ns('')^-Km('') 



(15) 



k = (— iy +l+1 / 2 (j + 1/2) and E n>K is the Dirac energy. For the point-charge nucleus, the 
sums £1 and £2, can be evaluated analytically, employing the method of generalized virial 
relations for the Dirac equation in a central field [7,8]. For k 7^ ±k' one can derive [7,8] 

|£i(k', n«)> ={[1 - (« - k') 2 ] [1 - (« + k') 2 ] + 4(qZ) 2 }- 1 

[!-(« + ^) 2 ] f + («' " ^ + r~'a z 



— (m e a x + E n)K io y )^j + 



4«Z[(k + At')m e - £„ iK ] 



K — K 



+ 2aZ[cr x r 1 + (ac + k')t 1 icr y ] \nk). 

and [11] 

|e 2 (K',n/«)) ={[4 - (« - k') 2 ][4 - (« + k') 2 ] + le(aZ) 2 }- 1 



(16) 



2«Z(k 2 - k' 2 )^ 



+ (« + k')[4 - (k - k') 2 ]5 + (8 + 8(aZ) 2 - 2(k - 



r 

8aZm e {K + k ) J — - + — 



2 Di 



(m e <Xr + E nK ia v 



j r D r k + k' Do y 
AaZm e D 1 2m e (n + k!) - 2E H)K D 2 16aZm 2 e (K + k') 



k 2 — k/ 2 -Dq 



K — K 



Do 



K — K 



where a y and ^ are the Pauli matrices, 



Do =[1 - (k + k') 2 ][1 - (« - ^) 2 ] + 4(«Z) 2 , 

£>i =2m e [8(aZ) 2 (K + k') 2 + (1 - (k + k') 2 )(8 + 8(aZ) 2 - 2(k - k') 2 ) 
+ 2£, ms (k + k')[8(«Z) 2 + (1 - (« + k') 2 )((k " «') 2 " 4)], 

D 2 =8aZm e [4 + 4(aZ) 2 - (k - k') 2 - (« + k') 2 (i - (k - k') 2 )] 
+ 12aZ J B niK ( K + «')[(« - k') 2 - 2]. 



|n«), (17) 

(18) 
(19) 
(20) 
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Further calculations of expressions (|1U[). 1)111) can easily be performed by using the re- 
current formulas for the expectation values A s = {nK\r s \n,K), B s = {nK\a z r s \nK,) , and 
C s = {nK\a x r s \nK) [7-9]. 

Finally, for the ns state, we obtain 

Q (ns,Z,I,») -/ me ^ f 1 ( 4E l-± + m * 

M ' I §(m e + E n -i) \ m e + E n -i m e - E n -\ 

| 3^ w ,_im e -2^_ ^ 2£ ra _i + m e -j ^ 
nij / 4m 2 , ' 



and 



Q Q (ns, Z, I, Q N ) 12Zm 2 (15 _ 16(aZ )2)( 45 + 4 ( aZ )2)( m 2 _ £ 2 _j 

10(-1728^ _ 1 - 3537m e #n,-i + 6237£ 2 _ 1 m 2 + 1809£?„ rl m] 
- 2457m 4 ) + m 2^ Z J 2 — (7405m^ - 736£g _ x + 10528m e ^_ 1 



e n,~ 1 

- 31033F 2 rl m^ + 2406^ _ x me - 3934^ rl mf + 20764E* ^m 2 ) 
+ T £ (6505m^ - 1120£^ > _ 1 m^ + 2240£ n rl mJ + 1280^ _ x 

m e — ^n,-l 

+ 7872m e Ei_ 1 - 1612^.^) + 400 ( TO g + 4^ )}. ( 22 ) 

It can be seen that formula (0) is a particular case (n = 1, F = 1, 7 = 1/2 (77 = 1)) of 
formula (jUJ). 

For small we can expand (|21j) and (1221) in the parameter aZ with the two lowest-order 
terms kept: 

Q> S ,Z,I,/i) = ^-^{{n 2 + 2) - {aZf{n - 1(1 - -L)) + 0((aZ) 4 )}g 1 . (23) 



r> /• C7 r /n \ t~Qtv f/61 4 355 2 , 
QQ[ns,Z, I,Q N ) = — ^-(Cgg " ~ + 4 ) 

^4(l028n 6 - 1830n 5 - 4322n 4 + 5325n 3 + 3957n 2 - 2808) + 0((aZf)\. 
540re^ ' ) 

(24) 



+ 



The main term in equation (|23|) coincides with formula (|4*|). 
For the state, a similar calculation yields 



1 ' bZ un { 6 V m e + E n ± (m e — E r _ 

}Qi- (25) 



5Z un L 6 V m e + E n> i (m e - -E^i) 2 
2E% 1 + 3m e E njl l \ 2E„. 1 — W-e 
m?_ m e — 1 / 4m,2 
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X 




12Zm2(15 - 16(aZ) 2 )(45 + 4(aZ) 2 )(m e - E n , x ) 
54 ( - + SllmeE^j + 6m^ n>1 - Ulm 3 e ) 



+ 



Qll2m e El A + W710E nA m 5 e + 211S2E^ 1 m 2 e - 998E% tl ml - 640£ , £ 1 ) 



+ 




For n = 2, F = 1, 1 = 1/2 formula (|25|) agrees with equation (jSJ). For small Z, we can 
expand expressions (|25jl and (|26|) in the parameter aZ' with the two lowest-order terms kept: 



3 Numerical results 

In Tables 1 and 2, we present the numerical results for Q^, Qq, and Qtotal = ^-Qn+Q/i+Qq 
(the total quadrupole moment of the atom) for a wide range of hydrogenlike ions. For ions 
with / = 1/2, the total quadrupole moment is completely determined by QM M . For the other 
ions, the role of the induced quadrupole moment is most important for low Z and decreases 
with Z increasing. 

We note also that, according to formulas (|23j) . (|27|). (|24[). (J2SJ, the induced quadrupole 
moment increases rapidly with n increasing. As a result, for highly exited states the total 
quadrupole moment of a hydrogenlike atom is mainly determined by the induced quadrupole 
moment. 
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Table 1: The numerical results for the QM M , the QMq, and Qtotai hi case of the Is state. The 
values of h//in and Qn are taken from [10]. 



Ion 


X H 


13Qt>+ 


17Q7 + 


43 Ca i9 + 


Z 


1 


6 


8 


20 


I 


1/2 


1/2 


5/2 


7/2 


fi/fi N 


2.79285 


0.702412(2) 


-1.8938(1) 


-1.3176 


Qn, barn 








-0.02578 


-0.049(5) 


Q li ,F = I + 1/2, barn 


0.7560 


0.03167 


-0.06401 


-0.01771 


Q Q ,F = I + 1/2, barn 








0.001072 


0.0008(1) 


Qtotai, F = I + 1/2, barn 


0.7560 


0.03167 


-0.08872 


-0.066(5) 


Qn,F = I - 1/2, barn 








0.02560 


0.009485 


Q Q ,F = I- 1/2, barn 








0.0008576 


0.0007(1) 


<2totai,-F = / - 1/2, barn 








0.00583 


-0.033(5) 



Ion 




207 Pb sl+ 




235TJ91+ 


Z 


54 


82 


83 


92 


I 


3/2 


1/2 


9/2 


7/2 




0.691862(4) 


0.592583(9) 


4.1106(2) 


-0.39(7) 1 


Qn, barn 


-0.120(12) 





-0.50(8) 


4.936(6) 


Q^F = 1 + 1/2, barn 


0.003281 


0.001697 


0.01158 


-0.0009(2) 


Q Q ,F = I +1/2, barn 


0.0007(1) 





0.0014(2) 


-0.01095(1) 


Qtotai, F — I + 1/2, barn 


-0.12(1) 


0.001697 


-0.49(8) 


4.924(6) 


Qn, F — I — 1/2, barn 


-0.0005469 





-0.007206 


0.0005(1) 


Q Q ,F = I- 1/2, barn 


0.00033(3) 





0.0013(2) 


-0.00978(1) 


Qtotai, F — I — 1/2, barn 


-0.06(1) 





-0.47(7) 


4.398(5) 



1 An average of the values given in [10]. 
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Table 2: The numerical results for the QM M , the QMq, and Qtotai in case of the 2py 2 state. The 
values of /x//xjv an d Qn are taken from [10]. 



Ion 


X H 


13Qb+ 


17Q7 + 


«Ca i9 + 


Z 


1 


6 


8 


20 


I 


1/2 


1/2 


5/2 


7/2 


fi/fi N 


2.79285 


0.702412(2) 


-1.8938(1) 


-1.3176 


Qn, barn 








-0.02578 


-0.049(5) 


Q^F = 7 + 1/2, barn 


-227170(100) 


-264.1 


300.0 


13.15 


Qq, F — I + 1/2, barn 








6.051 


0.73(8) 


Qtotai, F = 7 + 1/2, barn 


-227170(100) 


-264.1 


306.0 


13.83(8) 


Q li ,F = I -1/2, barn 








-120.0 


-7.044 


Q Q ,F = I -1/2, barn 








4.838 


0.65(8) 


Qtotai,-? 1 = J - 1/2, barn 








-115.2 


-6.44(8) 




Ion 


131 Xe ^+ 


207 Pb 81+ 


209gj«2+ 


235Tjyi+ 


Z 


54 


82 


83 


92 


I 


3/2 


1/2 


9/2 


7/2 


A*/ A* at 


0.691862(4) 


0.592583(9) 


4.1106(2) 


-0.39(7) 1 


Qat, barn 


-0.120(12) 





-0.50(8) 


4.936(6) 


Q M ,F = 7 + 1/2, barn 


-0.3096 


-0.06061 


-0.4015 


0.025(5) 


Q Q ,F = 7 + l/2, barn 


0.09(1) 





0.10(3) 


-0.723(1) 


Qtotai, F = I + 1/2, barn 


-0.34(1) 


-0.06061 


-0.8(1) 


4.24(1) 


Q il ,F = I -1/2, barn 


0.05160 





0.2598 


-0.013(5) 


Qq, F = I - 1/2, barn 


0.045(5) 





0.093(15) 


-0.645(1) 


Qtotaij-f = 1 ~ I/ 2 , barn 


0.04(1) 





-0.11(9) 


3.75(1) 



1 An average of the values given in [10]. 
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